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ABSTRACT 



One topic of current interest in thermoacoustic research is an annular prime 
mover [Lin et al., J. Acoust. Soc. Am, 100, 2846 (1996)]. The starting point for 
this research is an investigation of a constricted annular resonator. A literature 
search of the field resulted in surprisingly few references. The results of analytic, 
numerical, and experimental investigations are presented. Introducing a constriction 
into an annular resonator splits each longitudinal duct mode into two modes, one of a 
higher frequency with a pressure antinode at the constriction and one at a lower 
frequency with a velocity antinode near the constriction. The lower mode is more 
sensitive to changes in the length and porosity of the constriction than the higher mode. 
Overall agreement between measured and predicted mode shapes and resonance 
frequencies is veiy good. It was found that it is necessary to include end 
corrections at the constriction to get accurate agreement between measured and 
predicted results. 
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L 



INTRODUCTION 



The purpose of this thesis is to investigate acoustic resonances in a constricted, 
annular resonator. The motivation for this project is based on a related problem in the 
area of thermoacoustis - the annular prime mover. The intent is that the measurements 
reported in this thesis will serve as preliminary validation of models being developed for 
the prime mover project. 

The azimuthal modes of a uniform cross section annulus are degenerate. The 
introduction of a constriction splits the degeneracy resulting in two modes having different 
resonance frequencies and distributions of acoustic pressure (mode shape). The goal of 
this thesis is to both measure and predict the frequency splitting and mode shapes. 

A theoretical discussion will be covered in the Chapter II. Chapter in is a 
description of the experimental apparatus including the annular acoustic resonator and 
peripheral equipment. In addition, a detailed discussion of the data collection procedure 
will be covered. Chapter IV is an analysis of the data obtained and comparison with 
predictions. The thesis is summarized and important conclusions are addressed in Chapter 
V. 
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H. THEORY 



A. INTRODUCTION 

The purpose of this chapter is to outline the theoretical analysis of a driven, 
constricted annular resonator. The discussion begins with the solution of the wave 
equation in and the effective radius of a uniform cross section annular resonator. This 
solution is developed in cylindrical coordinates. However, because the purpose of 
this thesis is to investigate frequencies below those required to excite cross modes, we 
adopt a simplified plane model for the majority of analysis. 

The next step is to introduce a constriction. It is shown that the constriction 
results in frequency splitting of what were degenerated, orthogonal modes in the uniform 
cross section resonator. It is also shown that the degree to which a particular mode is 
excited depends upon the location of the driver relative to the constriction. Because of 
dissipation, the modes overlap in frequency space and, in general, two modes are 
simultaneously excited. When this situation occurs the resonance frequencies of the 
individual modes cannot be determined from peaks in the frequency response. 
However, pole-zeros analysis of the frequency response can yield accurate values for the 
resonance frequencies. The chapter concludes with a discussion of pole-zero analysis of 
two coupled modes. 

B. SOLUTION FOR A UNIFORM ANNULAR RESONATOR 

The analysis begins with the general solution of the acoustic wave equation 
in cylindrical coordinates 




( 2 . 1 ) 



1 




1 d 2 d 2 



( 2 . 2 ) 



where 




r 
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Assuming harmonic time dependence, the wave equation becomes the Helmholtz 
equation 

V 2 P = -k 2 ? where k 2 =W . (2.3) 

The general solution of Eq. (2.3) is assumed to be expressed as the product of three 
terms, each a function of only one variable. 




Figure 2-1. Geometry for the annular resonator. 



P (r,6,z) = R(r)Q>(0)Z(z) 

Substituting Eq. (2.4) into (2.3) and dividing by R 0 Z , gives 



J_J__£_r dR \ 
R r dr v dr J 



1 1 d 2 6 



0 d6 



1 dj_ Z 
Z dz 2 



= - k 



Rearranging to isolate the z - dependent terms gives 



1 1 d 



R r dr 



dR 

r 



\ 



dr 



1 1 d 2 0 



+ k 2 = - 



1 d 2 Z 



r 0 d 0 1 Z 

Following the standard separation of variables solution technique , 

1 d 2 Z 



dz 



let 



= - k 



Z dz 2 

The general solution to this equation is 

Z - A. cos ( k 2 z ) + B z sin(Ar z z) 

Equation (2.6) now becomes after multiplying through by 
rearranging 



(2.4) 



(2.5) 



( 2 . 6 ) 



(2.7) 

( 2 . 8 ) 

and 
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r_ £_ 

R dr 



r Sr) + (* 2 - */>’ = 



1 <? 2 0 
0 dQ 2 ‘ 



Let 



1 2 0 



= - m 



0 dG 2 

The general solution to this equation is conveniently expressed in the form 

0(0) = cos (m 0 + ^ ) 

Equation (2.9) is now reduced to 



R dr 






+ (k 2 r 2 - m 2 ) = 0, 



where 



k 2 - k 2 - k . 



(2.9) 

( 2 . 10 ) 

( 2 . 11 ) 

( 2 . 12 ) 

(2.13) 



The general solution to this equation is 

R(r)= A m J M (k,r) + B m r m (k,r), 



(2.14) 



where J m (k r r) and Y m {k r r ) are the cylindrical Bessel and Neumann functions. The 

Neumann function must be retained because the origin is excluded in our annular 
resonator. 



All boundaries are assumed to be rigid, so the normal component of acoustic 
velocity must be zero at the boundaries. Applying rigid boundary conditions at z = 0 
and z - L 2 implies that B_. = 0 and 



, / =0,1,2... (2.15) 

The function 0(0) must be periodical, in other words, 0(0) = 0(0 + 2k) . This 
condition is satisfied if m is an integer, m = 0,1,2 ••••. 

Application of rigid boundary condition at r = a and r = b implies that 

A m r m (M) + B m r m (M) = o 

A m J n (k r b) + B m Y m (M)= 0. (2.16) 

Expressed in matrix form these equations become 
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(2.17) 



The values of k r a 
first matrix is zero , 



J'iKb) 



f°l 

vOy 



that satisfy this equation are those for which the determinant of the 



i.e. 



J m (k r a)Y m (k r b) - J' m (k r b)Y' m (k r a) = 0 . (2.18) 

Let the values of k r a that satisfy Eq. (2.18) be called a mn . With this definition, the 
first of Eqs. (2.16) becomes 



and Eq. (2.14) becomes 



B 



= -A 



Y m (a mn ) ’ 



R(r)= A m [J m (k r r)~ 



y 

Y m («„») ' 



(Kr)l 



(2.19) 



( 2 . 20 ) 



For the purposes of this thesis, we are interested in frequencies low compared 
to those required to excite cross modes. To eliminate pressure variations in the z 
direction, we require that / = 0 in Eq. (2.15). In this case Z(z) = 1 , and k ? = k 2 . 
We are also interested in modes for which the circumference of the resonator is 
approximately one wavelength long, i.e. m = 1 in Eq. (2.11). Under these conditions, 
the solution to the Helmholtz equation is 

P (r,$) = A [./,(*,/•) - r'!~ , '4 |, i( t . f )] |:os ^ + (2-21) 



The dimensions for our particular resonator are a = 10.00 cm and b = 15.30 cm. The 
lowest value of a mn is then a 10 = 0.796. Recall that for 1 = 0, 

k l = k 2 = ( --— ) 2 . Also, k r = . Combining these expressions for k r 



implies that 



_ & ip 

X a 

Let us define the effective circumference of the resonator L eff 



( 2 . 22 ) 

m X , or for m = 1 



L eff = X . Combining this definition with the Eq. (2.22) and rearranging yields 



6 



2 na 



(2.23) 



'ff 



a 



10 



Alternatively, we can define an effective radius R eff = . 



Using Eq. (2.23) , we 



a 



have R eff = — — . For a I0 = 0.796 , R eff = 12.56 cm. 



a 
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C. PLANE WAVE THEORY 



The parameter R eff is defined in such a way that 



R eff ~ 



‘eff 



k 



In 2 n 

The angular coordinate 0 can be related to the arc length x measured at R eff 



(2.24) 

through 



6 = — = kx . (2.25) 

Reff 

With this definition , the 0 - dependence of the acoustic pressure becomes 

cos (kx + <j>). (2.26) 

This function can be expressed in terms of complex exponentials as 

cos (kx + </>) = 2[e' ( ** + * ) + e - JVcx+ * ) ] , (2.27) 

Substituting this result into Eq. (2.21) gives 

P (r,e,t)= R(r)[e K ‘ 0, - kx -* ) + e ><•'+**♦*>] ( 2 . 2 8) 

The interpretation of Eq. (2.28) is that the pressure is represented by two counter 
propagating waves that behave essentially like plane waves. This analogy will be 
used to analyze the nature of the sound field in the annular resonator. 



The plane-wave-like solution will be valid for frequencies below those required 
to excite the first cross mode. The frequencies corresponding to normal modes in 
the r - z plane are given by 



7 



(2.29) 



fc = 






Table 2-1 lists a few values of a, 



m\n 


0 


1 


2 


1 


0.796 


6.025 


11.903 


2 


1.589 


6.119 


11.987 


3 


2.378 


6.465 


12.124 



Table 2-1. A few values of a, 



In our annulus a =0.10 m, so it is seen that 



a. 



is approximately 60 for m = 



K 



1. Similarly, L z is 0.05 m, so — is equal to 62.8, or nearly equal to 



a 



1 If 



a 



Therefore the frequencies of the first cross mode in either the z or r direction are 
approximately equal and have value in the vicinity of 3.3 - 3.6 kHz. Our primary 



interest is in frequencies corresponding to = A , or / = — — = 440 Hz. 






1. Uniform Cross Section Annular Resonator 
a. Undriven case 

As explained above, for low frequencies the annular resonator can be 
modeled as a straight pipe of length L ^ with periodic boundary conditions at the 

two end as shown in Fig. 2-2. Assume that there are two traveling waves in pipe. 
One is propagating in the positive x-direction and represented by 

P + = Ae Ka, - kx \ (2.30) 

The other propagates in the negative x- direction and given by 

P_ = Be J(a>,+kx) . (2.31) 
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The total pressure in the pipe is a sum of these two pressures and given by 

P ( = P + + P_ = Ae + Be * at+kx) . (2.32) 

These subscripts get used later in a different way use p and n for positive and negative 
instead of 1 and 2. 




Figure 2-2. Annular resonator and equivalent straight resonator. 



Also, the volume velocity U can be obtained from knowledge of the pressure and the 



pc 

pipe’s acoustic impedance, Z = — — , where S is cross section area of pipe, 

J 

Sk 



U t = U + - U. = [Ae J ' (<B, - te) - B j(a,t+te) ] . 



Periodic boundary conditions are 



and 



ap i 



P = P 

r x=0 r x=L 



Ux. 0 = U X ._ L 



where L = L eff . Substituting Eq. (2.32) into (2.34) gives 

e* - 1 



(2.33) 

(2.34) 

(2.35) 

(2.36) 



Substituting this expression into Eq. (2.33) and applying Eq. (2.35) yields 

cos kL = 1 . (2.37) 



This implies that 



kL = 2nm 



( m = 0 , 1 , 2 ), 



(2.38) 
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or 




(2.39) 



It should be pointed out that L eff has slightly different values for different values of 



m. 



b. Driven case 



In addition to knowing the resonance frequencies of the annular 
resonator, we are also interested in the mode shapes, i.e. the acoustic pressure and 



a nonuniformity in cross section must be introduced to fix the mode shapes. In the 
undriven uniform cross section case, any superposition of two degenerate 
orthogonal modes ( one proportional to cos mO , the other to sin md ) satisfies the 
boundary conditions. In order to find the mode shape, it is necessary to install a 
driver in resonator. Assume that a branch driver is located at x = 0 as shown on Fig 
2-3. The volume velocity injected by the driver is given by 



velocity as function of angular position. There is no unique orientation of the modes 
in an undriven, uniform cross section annulus. Either the system must be driven, or 




(2.40) 



x=0 



x=L 




Figure 2-3. Driver installed straight resonator. 



The boundary condition for pressure is given by 

P = P 

X x=0 r x=L 

which is same as that for the undriven system. This conditions leads to 



(2.41) 




(2.42) 
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The boundary condition for volume velocity now includes an additional term and is given 
by 



which yields the equation 



U 



x=0 



~U^+U D , 



(2.43) 



A(l-e-*)-B(l + ^)=^ I) . (2.44) 



Equations. (2.42) and (2.44) can be expressed in matrix form as 



'l - e~ jkL 


1 - e JkL ' 


"A" 




0 


1 

1 

1 


- 1 + e jkL 


B 


— 


Po C 

S U D 



(2.45) 



This equation has the form MC = U . The value of the coefficients A , B can be 



determined from C = M l U . Once A and B are known, the pressure and volume 
velocity are found from Eqs. (2.32) and (2.33). 



2. Constricted Annular Resonator. 

a. Frequency vs area ratio (no driver) 

The analysis of the constricted resonator is more complicated due to its 
non uniform cross sectional area. The geometry is shown in Fig. 2-4. 



si 

1-* w 


S2 


R W 






x=0 x= / x=L 

Figure 2-4. Constricted resonator with no driver. 



Calculation of the resonance frequency of a constricted resonator is similar to that 
of the uniform cross section resonator. In order to solve this case, the basic equations 



11 



of pressure and volume velocity in each section have same form as the uniform cross 
section resonator . 

P, = Ae i(at - k ' x) + Be * at * k ' x) , (2.46) 

P 2 = Ce Ka,t - k ' x) + De Jia, + kiX) . (2.47) 

At the junction of two sections, pressure and volume velocity are continuous which gives 
the following equations 



P,(i = /) = P J (x = /) > 


(2.48) 


P,(* = 0) = P,(x = L) , 


(2.49) 


II 

b" 

n 

II 

2 1 


(2.50) 


</,(*« o) = £/,(* = £)• 


(2.51) 



These equations can be combined into a 4 x 4 matrix in the following form. 



e ~ jkl ‘ 


e jk '‘ 


1 

1 


- e jk ' 1 ‘ 




1 


1 


_ e -jk 2 L 


_ e M 


' A ' 


e 


_ e M 


~ ^2^2 c ~jk,l 


$2 k 2 jk,l 


B 






Stk t 


S ] k l 


C 


1 


- 1 


~ ^2^2 c -ik,L 


$2 k 2 ik,L 


_D 






S \ k \ 


S x k x \ 





= 0 • (2.52) 



This matrix has the form of M C = 0 . The resonance frequencies are those 
frequencies that satisfy the equation det[M] = 0. 

Unlike the uniform cross section case, the introduction of a constriction 
splits the degeneracy of the modes. The two modes, called the high frequency and 
low frequency modes, correspond to different mode shapes. ( In what follows, the 
word “frequency” be will frequency dropped and the modes referred to, simply, as 
“high” and “low”.) The high mode has a pressure antinode centered on the constriction 
whereas the low mode has a velocity antinode centered on the constriction. The 
difference in frequency between the high and low mode increases with increasing 
length or decreasing cross sectional area of the constriction. 
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The dependence of the frequency splitting is shown in Figs. 2.5 a-c. 
The plots show resonance frequency as a function of area ratio for three different 
constriction lengths. The resonator contains air at 20°C and atmosphere pressure. Area 
ratio is defined as the ratio of the cross sectional (open) area in the constriction to 
that of the resonator. In the limits of zero area ratio, i.e. when the constriction is 
completely blocked, the high mode frequency corresponds to X = (L eff - L comtriction ) 

and the low mode frequency corresponds to X = 2(L <ar - L constriction ) . These limits are 
appropriate for the two lowest modes of a rigid, rigid resonator of length - L constriction . 



calculated resonance freq vs area ratio 




Figure 2-5a. Resonance frequencies with a long constriction 



t \ 

16 J 
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calculated resonance freq vs area ratio 




Figure 2-5b. Resonance frequencies with a medium length constriction 



132 ; 



calculated resonance freq vs area ratio 




Figure 2-5c. Resonance frequencies with a short constriction 
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b. Mode shapes of a constricted resonator 

Calculation of the pressure distribution in a driven, constricted resonator 
is also similar to that for the driven uniform cross section case. For the driven, constricted 
resonator, the resonator is divided into three sections, as shown in Fig. 2-6. The driver 
is located at x = d. 



Si | S2 | S3 

< ►}< ► 




x=d 



Figure 2-6. Constricted resonator with driver. 

This system gives rise to the following set of equations for pressure and volume velocity in 
the three regions: 



P, = Ae J(a,l ~ k ' x) + Be i{a,+k ' x) , 


(2.53) 


P 2 = Ce K *‘- klX) + De j(at+k ' x) , 


(2.54) 


P 3 = Ee J( - a, - ktX) + Fe KM+k2X) , 


(2.55) 


U 1 ^ 1 ^ 1 Ac l(e>t ~ k,x) 

' c op a cop 0 


(2.56) 


U ^2^2 QgKaH-k-yx') $2^2 j(a>t+k 7 x) 

' cop 0 ^ (op 0 


(2.57) 


e l C Ir 

JJ °2 K 2 ]? c not-k,x) °2 k 2 p c j{<Dt+k. 1 x) 

3 (Op o ^ cop Q 


(2.58) 


U D = U D e Ja * ■ 


(2.59) 



The three different sections also give rise to six boundary conditions: 
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P,(x = 0) = P s (x = I), 


(2.60) 


U,(x = 0) = U 3 (x= L), 


(2.61) 


II 

aT 

ii 

ii 

aT 


(2.62) 


U,(x = t) = U 2 (x = l), 


(2.63) 


P ; (x = <#) = Pj(x = d), 


(2.64) 


U 2 (x = d) + u D = U,(x = d). 


(2.65) 



Application of these boundary conditions to Eqs. (2.53) - (2.59) yields a set of equations 
which can be combined into a matrix equation of the form MC - U . 



1 


1 


0 


0 


- e~ ]1ClL 


1 

< 

CD 

1 


* A ^ 




0 


1 


-1 


0 


0 


— c k 


^2^2 jk 2 L 


A 

n 




0 










S : k } 


SA 


B 






e~ ik ' 1 


e*' 1 


-e~ Jhl 


-e*> 


0 


0 


c 




0 


e~** 


- e*' 1 


-SjK -jm 




0 


0 


D 


— 


0 








S } k, 






p 




0 


0 


0 




e Jhd 


-e~ M 


-e** 


JZf 




~Po C 


0 


0 


e~** 


-e** 


-e** 


e** 


F 




Ud 



Solving the inverse form of Eq. (2.66) gives the value of the coefficients A through F. 
Substitution of these coefficients into Eq. (2.53)-(2.55) gives the mode shapes. 

(1) Mode Selection and Driver Location. Because of 
dissipation, the quality factor (Q) of a particular mode is finite. This means that, in 
general, when the resonator is driven at a frequency in the vicinity of the resonance 
of either the high or low mode, both modes will be excited. In order to measure 
the mode shape of a single mode, it is necessary to excite only that mode as nearly as 
possible. This can be accomplished to an acceptable degree by proper location of 
the driver. If the driver is positioned at a pressure antinode of a particular mode 
and driven at the resonance frequency of that mode, then that mode is the 
predominant one excited. 

Mode selection is illustrated the next sequence of figures which show the 
calculated frequency response and mode shapes for the resonator with a constriction 
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length of 



and area ratio 0.7. For Fig. 2-7, the driver is located 180° from the 



L eff 

32 



center of the constriction which should correspond to a pressure antinode for the 
high mode. Figure 2-7a shows the frequency response as indicated by computed 
magnitude of acoustic pressure measured 45° from constriction. This figure shows that 
there is only a single peak in the frequency response. The peak frequency, 439.0 Hz, 



corresponds to the resonance frequency of the high mode for a 




long constriction 



with a 0.7 area ratio, as shown in Fig. 2-5b. Figure 2-7b shows the calculated mode 
shape when the resonator is driven at 439.0 Hz. It is seen that there is indeed a pressure 
anitnode at both the constriction and driver location. It is also evident that the standing 
wave ratio is very high, another indication of single-mode excitation. 



Frequency Response 




Figure 2-7a. Frequency response of the constricted resonator when driver is located 1 80° 
from the constriction and driven at the high mode resonance frequency for the 




length, 0.7 area ratio constriction. 
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Mode Shape for Length L/32, Area Ratio 0.7, and f=439.0 




Figure 2- 7b. Mode shape when the resonator is driven at the peak frequency as 
determined from Fig. 2-7a. 

Figure 2-8a shows the frequency response when the driver is moved 
to a location 90° from the constriction. This driver location should be near a 
pressure node for the high mode and pressure antinode for the low mode. The 
frequency response again shows only one peak. Comparison with Fig. 2-5b indicates 
that this mode is the low mode. Figure 2-8b shows the mode shape. It is seen that 
now there is a pressure node at the constriction. Again, the standing wave ratio is 
high indicating single-mode excitation. 
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Figure 2-8a. Frequency response of the constricted resonator when driver is located 90° 

A# 

from the constriction and driven at the low mode resonance frequency for the 
length, 0.7 area ratio constriction. 



Frequency Response 




Mode Shape for Length L/32, Area Ratio 0.7, and f=427.2 




Figure 2-8b. Mode shape when the resonator is driven at the peak frequency as 
determined from Fig. 2-8a. 
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In the cases of the driver located at 180 degrees or 90 degrees from 
the constriction, the driver excites only one mode, either high or low. If the driver is 
placed at a position other than 1 80 degree or 90 degree, two modes will be excited 
simultaneously. Figure 2-9 shows the frequency response when the driver is located 45 

degrees from the constriction. We have used the 0.1 area ratio constriction to 

make the amount of frequency splitting larger. 



Frequency Response 




Frequency (Hz) 

Figure 2-9a. Frequency response when driver is located 45° from the ^/j^, 0.1 area 
ratio constriction. 

We can see from Fig. 2-9a that both the high and low modes are excited. In this case, 
they are well separated, in other words, there is no significant overlap. When there 
is little overlap, it is easy to pick up the resonance frequency of each mode. A constriction 

of small area ratio (0.1) and long size (^/^, 22.5 degree), as in Fig. 2-9a, separates the 

two modes well enough for there to be little overlap. Figure 2.9b shows the mode 
shape when the resonator is driver at the lower peak frequency. This figure shows 
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that even when the driver is at a location which would allow both modes to be 
excited, when the modes are well separated, single mode excitation is still possible. 



Mode Shape for Length L/16, Area Ratio 0.1, and f=320.0 




Figure 2-9b. Mode shape when the resonator is driven at the lower of the two 
peak frequencies from Fig. 2-9a. 



Figure 2- 10a shows the frequency response for the 




length, 0.7 area ratio constriction used in earlier figures. In contrast to Fig. 2-9a, the 
modes show significant overlap. This overlap causes two problems. First, it means 
that a driver located at 45° will excite both modes to a significant degree, even 
when driver at one of the peak frequencies. Second, the peak frequencies no longer 
correspond to the resonance frequencies. Overlapping is most significant for large area 
ratio, short constrictions. 

Figure 2- 10b shows the mode shape when the resonator is driven 
at the lower peak frequency in Fig. 2- 10a. It is seen that the standing ratio is low in 
this case, indicating excitation of more than one mode. The next section addresses a 
technique to determine resonance frequencies in the case of significant overlap. 
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Frequency Response 




Figure 2- 10a. Overlapped wave form of frequency response with driver 45 degrees from a 
constriction which has a 0.7 area ratio and (5.625 degree) length. 



Mode Shape for Length L/32, Area Ratio 0.7, and f=428.5 




Figure 2- 10b. Mode shape when the resonator is excited at the lower peak frequency 
in Fig2-10a. 
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D. 



POLE -ZERO ANALYSIS 



For understanding the acoustics of an annular resonator, it is important to find 
accurate values of the resonance frequencies. These resonance frequencies are 
independent of the driver location but are dependent upon what kind of constriction is 
being used - the length and area ratio. As we have seen, the driver may excite one or 
two modes simultaneously depending upon its excitation frequency and location. When 
only one mode is excited, it is simple to find resonance frequency by finding the 
frequency corresponding to the peak in the frequency response . But for the case when 
the driver excites two modes simultaneously, the peak frequency cannot be used to 
determine the resonance frequency. This situation was discussed in the previous section. 
Pole-zero analysis can be used to determine the resonance frequencies in the case of 
significant modal overlap. 



1. Analysis 
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For a driven simple harmonic oscillator, the complex speed is represented by [Ref 
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The expression for H can be written as 
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by introducing the resonance frequency co a and quality factor, Q = ( ° 0 /r ■ 

/ ftt 

By analogy, we assume the pressure in a resonator driven near resonance, at 
one point in space, has the form 
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where Q is now the quality factor for the resonator. Factoring the denominator of Eq. 
(2.69) results in 
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This equation can be simplified if the denominator is written as [Ref 2] 
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where the poles p x and p 2 are. 
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where co l = co 0 , for the single mode case. The poles are written in the form (co + jp) 
because the Hewlett Packard signal analyzer used for experimental measurement outputs 
the poles in that form. Equation (2.74) can be simplified into a form given by 



p--a± jb where a = 
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For two overlapping modes, the complex pressure amplitude is the sum of the individual 
amplitudes 



= Pi + P; 



(276) 



and P Iot has the form 

p = A, +Jp 1 Xo) +jp ± l+ A 2 (a> + JPi)(<a ±JPi) 
,ot (a+jPi ) (a + jP 4 ) 



(2.77) 



So, even though the two modes are overlapped and the amplitude of the individual 
responses are added, the poles still show separable frequencies. 

The numerator of Eq. (2.77) can be factored in the same way as the denominator 
of Eq. (2.69). In this case the numerator is 
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From Eq. (2.78) we can find zeros. 

The resonance frequency and quality factor can be obtained from the poles using 
Eq. (2.75). The resonance frequency is then 

a), = a/<7 2 + Z> 2 (2.81) 
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and the quality factor is 



co, 

Q . = Ya ' 

2a 



(2.82) 



co. 



If Q is large, an approximation can be used to simplify Eq. (2.81). The term of — r - 

2Q\ 



can be neglected since 



\ 2 
co. 



U 0, 



is small compared to co , 2 . So, b = co , . 



For data analysis, we fit the measured frequency response using pole-zero 
analysis. The fit parameters give the resonance frequency and Q of the individual modes. 



THERMO- VISCOUS DAMPING 



Up to this point, losses have been ignored. However any real gas has viscosity 
and will exchange heat with the walls of the resonator. The boundary condition at 
the walls of the resonator is that the normal velocity is zero. When viscosity is included 
the correct boundary condition has both the normal and tangential velocity going 
to zero. This creates a viscous boundary layer which contributes viscous losses to the 
propagation of sound [Ref. 3,4], Similarly, the gas near the walls will exchange 
heat with the walls. This thermal boundary layer also contributes losses to the sound 
propagation [Ref. 3,4]. While these viscous and thermal effects are complicated, if the 
boundary layer thicknesses are small compared to the cross sectional dimensions of the 
resonator, the effects can be included by a modification of the wave number k. 

The frequency dependent viscous and thermal boundary layer thickness 8 V and 



8 k are given by 



and 




(2.83) 

(2.84) 
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where p is the viscosity, K is the thermal conductivity, c p is the specific heat and p 
is the density of the gas. The wave number k is then given by [Ref. 5] 



where r h is the hydraulic radius, which is the ratio of the cross sectional area to the 



and in the main channel, and that k is now complex. The imaginary part of k is related 
to the thermo-viscous damping. 

F. END EFFECTS 

In the earlier analysis, the boundary conditions used at the ends of the constriction 
were continuity of pressure and volume velocity. Using only plane wave analysis, the 
pressure and particle velocity are constant across the cross section. However, at the 
constriction end, the velocity and pressure have a spatial dependence and thus the plane 
wave theory is not as accurate near the ends as it is away from them. 

The spatial variation in the pressure can be described using non-propagating 
higher order modes. These non-propagating modes store energy near the constriction 
ends. Their effect on the plane wave propagation can be modeled as a lumped impedance 
element. [Ref. 6], The lumped impedance can be obtained from a variety of methods. 
A first approximation involves modeling the impedance as an added mass much like 
the end correction for radiation at the end of a tube. This added mass can be computed 
using conformal mapping. [Ref. 6], A more accurate method involves including higher 
order modes near the constriction end. [Ref. 7]. While computationally more intensive 
than the conformal mapping solution, the higher order mode solution is more accurate. 

After computing the lumped impedance Z (<y) , the boundary condition for the 
pressure becomes (the volume velocity boundary condition is unchanged.) 




(2.85) 



perimeter. It should be noted that since the hydraulic radius is different in the constriction 



P,(/) = P 2 (/) + Z (o>)U 2 (l) 



( 2 . 86 ) 
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where P, is the pressure in unconstricted region and P 2 is the pressure in constriction. 
U 2 is the velocity in the constriction. The constriction is at x = l. 
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m. EXPERIMENTAL APPARATUS AND PROCEDURE 



A. INTRODUCTION 

A description of the experimental apparatus and the procedures for data 
acquisition are presented in this chapter. The first section will include a description of the 
annular resonator, constrictions, the driver, the microphone and electronic 
instrumentation. The second section discusses the data acquisition system. The third 
section is a detailed discussion of the experimental procedure. 

B. EXPERIMENTAL APPARATUS 

1. Resonator and Driver 

Figure 3-1 is an illustration of the annular resonator and acoustic driver. The sides 
and bottom of resonator are made of aluminum. The top of resonator is made of 
plexiglas to allow visual sighting of the relative positioning of the constriction and the 
driver. More details about each part will follow. 

The dimensions of the resonator were chosen to give longitudinal resonance 
frequencies in the 300 ~ 500 Hz range and to be well below the cutoff frequency 
for cross modes. The resonator consists of two main parts as shown on Fig. 3-1. 
One piece is an aluminum ring with 15.30±0.02cm inner radius and lcm thick wall. 
This ring formed the outer side wall of the resonator. The other piece is a solid 
aluminum cylinder which is concentric with the outer ring. The inner radius of 
resonator is 10.0+0.01 cm. The height of resonator is 5.0+0.01 cm. The resonator 
is capped on the open side with lcm thick plexiglas plate with a radius of 16.0 cm. 
The bottom plate of the resonator is made of aluminum and its thickness is 1.0 cm. 
Sixteen microphones are installed in the cap to measure the spatial acoustic pressure 
distribution inside resonator. The installation of microphones will be discussed later. 
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The cap is bolted with 10.7 cm long threaded rod in the center and 8 clamp bars as 
shown on Fig. 3-2. To ensure a tight seal (for a high quality factor), vacuum grease 
was applied to the surface of the aluminum ring and center piece. 




Figure 3-1. Illustration of resonator and driver. 



The acoustic signal was generated by a 40 W , 16 ohm external driver. The 
driver was connected to the resonator by a 3.5 cm long piece of plastic hose with 0.4 
cm diameter. The driver was excited by the output of a Hewlett Packard 3 562 A dynamic 
signal analyzer sent through a Techron 5507 power amplifier. The signal analyzer 
was used to determine the frequency response of the resonator. Typically, the 
resonator was driven with random noise withal Vrms source level. To measure the 
mode shape, a Hewlett Packard 33120A function generator was used to excite the 
resonator at the resonance frequency determined from the built-in pole-zero analysis 
of the frequency response using the signal analyzer. 
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Figure 3-2. Illustration of clamp bar, threaded rod and nut. 



2. Microphones 

For the experiment, 16 Panasonic WM60-AY microphones were used to 
measure the pressure distribution inside resonator. The dimensions of the microphones 
are 6 mm in height and 5mm in diameter. The sensitivity is -42 dB ± 3dB re IPa. A 
Hewlett Packard 623 7B triple output power supplier provided 9 V DC to each 
microphone. 



a. Installation of microphones 

The microphones were installed as shown on Fig. 3-3. The first 
intention was to use only one microphone and rotate the cap for each different position 
to measure the mode shape. However, the rotation of cap broke the sealing between the 
plexiglas cap and aluminum bodies of the resonator making the measurements very time 
consuming and nonrepeatable. 
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mic 1 




Figure 3-3. Installation of 16 microphones on the plexiglas cap. 

We decided to install 16 microphones in order to protect sealing and reduce the time 
required for data acquisition. For the precise measurement in constricted area, 6 
microphones are located in first 30 degree region where a constriction is located. The 
remaining 10 microphones are located 30 degrees apart. The microphones were installed 
in the center of the annular channel. To install the microphones, 16 holes were drilled 
part way through the plexiglas cap. These holes were slightly larger than the diameter of 
the microphone. A smaller hole was drilled concentric with each of the larger holes. This 
smaller hole penetrated the cap. Therefore the hole size presented to the inside of the 
resonator was much smaller than the microphone diameter. The small hole diameter is 
0.15 cm. The microphones were sealed into the holes with silicon glue. The error in 
the microphone position is estimated to be less than ± 0.25°. 
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Built in speaker 
Chamber 




Figure 3-4 . Chamber used for microphone calibration 

b. Calibration of microphones 

Because 16 microphones are used for sensing the pressure distribution, the 
first step was to calibrate all the microphones over the 300 - 500 Hz frequency range 
interest. For the calibration, a custom made chamber with a built in speaker was 
used. As shown on Fig. 3-4, the diameter of the chamber is 3.0 cm and the height is 2.50 
cm. At first we tried to measure the pressure for each microphone at the same 
position relative to the driver in an otherwise empty resonator. However, the phase 
shift caused by the resonators’ resonance made calibration very sensitive to small 
errors in position. Instead, the top cap has removed and the calibration chamber 
was positioned directly over each microphone. The frequency response was measured 
for each microphone over a 300 ~ 500 Hz range. A typical frequency response is 
shown in Fig. 3-5. This sweeping and data storage were done by Hewlett Packard 
3 5665 A dynamic signal analyzer. For the calibration, the signal analyzer was set to swept 
sine mode with a settling time of 20 cycles and integration time of 50 cycles. The 
chamber was placed onto the cap where each microphone is located and the chamber 
was sealed on the cap with vacuum grease. 
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Y:474.0161 m 



OX:353 Hz 




Figure 3-5. Frequency response of Mic. #1 covered by the calibration chamber. 

Microphone #1 was used as the reference for comparison calibrations of the 
other 15 microphones. Each microphone’s output was divided by that of 
microphone #1 to create a correction factor as shown on Fig. 3-6. When measuring 
the mode shape, the output from each microphone was divided by the appropriate 
correction factor to obtain the corrected amplitude at a particular location and 
driving frequency. For example, consider the case when the measured amplitude from 
microphone # 12 is 50 mV at 455 Hz. The correction factor is 0.739 from the Fig. 
3-6. Dividing 50 mV by 0.739 gives 67.66 mV as the corrected measurement. 
Generally, the correction factors were between 0.7 and 1.2. 

O X:455 Hz Y:738.7183m 




300 Hz 500 Hz 



Figure 3-6. Plot of correction factor for Mic. #12 referenced to Mic. #1. The correction 
factor is 0.739 at 455 Hz. 
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3. 



Constriction 



To explore a wide range of parameters, various sizes of constriction are necessary. 
When considering the material selection for this constriction, a couple of factors came 
to mind. The material had to be rigid and have good machining propertied for ease 
of fabrication. These two factors led to the selection of PVC (Polyvinylchloride). 
Figure 3-7 shows the design along with the nominal and measured dimensions for the 
various constrictions used. For the experiment, 3 different angular widths were used. 
The width of constriction is given by the 0 shown in Fig. 3-7. The three widths 

are 22.5° , 11.25° and 5.625°. The 22.5° constriction corresponds to of the 

effective circumference of the resonator. 1 1 .25° corresponds to of L eff and 

5.625° corresponds to — of . The area ratio is the ratio of open area in the 

64 

constriction (which is the shaded black square in Fig. 3-7) to cross sectional area of 
resonator. In Fig 3-7, Si stands for the open area and S2 stands for the open cross 
sectional area of the constriction. The table in Fig. 3-7 gives the geometrical parameters 
of 32 different constrictions. 
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Figure 3-7a. Nominal geometry of the constrictions 



unit: cm 
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Figure 3-7b. Measured geometry of the constrictions. 



4. Electronic Instrumentation 



Figure 3-8 is a block diagram of the electronic instrumentation used during 
the measurements. A Hewlett Packard 3562A Dynamic signal analyzer is used to 
determine the frequency response in the band from 300 Hz to 500 Hz. The signal 
from this analyzer is amplified by a Techron 5507 power supply amplifier and used as 
the input signal to the driver. The pole-zero analysis function built into the signal 
analyzer is used to extract the resonance frequencies from the measured frequency 
response. Once the resonance frequency is found, a Hewlett Packard 33120A 
Function/ Arbitrary Waveform Generator is used to generate a 1 Vpp sinusoidal 
signal at the obtained resonance frequency. The signal is also amplified by the 
same Techron amplifier and used as the input signal to the driver. Either the signal 
from the analyzer or the function generator signal is monitored on a Tektronics 2336 
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Oscilloscope. The outputs from the microphones are monitored on the same oscilloscope. 
A custom made microphone switch board selects each microphone to measure the 
pressure amplitude at different positions. The microphone output amplitude is measured 
using a Standford Research 530 Lock In amplifier. In order to measure the 
temperature of the resonator, an Omega 450 AET thermocouple thermometer type E is 
used. The thermocouple was calibrated to ±0.1 °C using a distilled water ice bath. The 
thermocouple is positioned on the surface of the outer aluminum ring of resonator. 
The temperature inside resonator is measured frequently and be recorded whenever there 
is variation. 
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Figure 3-8. Block diagram of instrumentation setup. 





c. 



EXPERIMENTAL PROCEDURE 



1. Resonator Setup 

Prior to any data acquisition, the resonator is assembled with the constriction in 
the desired position relative to the driver. For the experiment, three major positions are 
used. They are 45°, 90° and 180° from driver. The position is measured from the 
driver to the center of constriction. The top cap is oriented so that 0° is set at the 
left edge of constriction. The 0° setting is performed by rotation of the cap by hand. 
For the error measurement of 0° setting, three people conducted rotation of the 
cap. Each person tried three times and we found that the error by rotation was less 
than ±0.5°. So, with the mechanical error ±0.25° in drilling the microphone hole, the 
total error of microphone positioning is less than ±0.75°. With the constriction 
positioned, the proper amount of vacuum grease is put on the open surface of resonator 
and the top cap is secured with the clamp bars. 

2. Determining Resonator Characteristics 

The Hewlett Packard 3 562 A Dynamic Signal Analyzer is used to characterize 
the resonator. The swept sine mode of the analyzer provides a sine wave signal to 
the system at a given range of frequencies. The frequency response of the system is 
used to determine the resonance frequency and quality factor of several modes of the 
resonator. When only one mode is excited, the peak pressure amplitude frequency is 
same as the resonance frequency given by pole-zero analysis. When two modes are 
excited, pole-zero analysis is used to determine the resonance frequency. The signal 
analyzer gives poles in the form a + jb . The resonance frequency is simply the 

value of b . The quality factor can be also obtained by from that form. This 

method has been explained in the theory chapter. These measurements are performed 
for each constriction type and each relative position of constriction and driver. 
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3. 



Data Collection 



Data collection is performed in sets of the positions and different sizes of 
constriction. Before data collection, the initial ambient temperature and noise level 
is checked. The output of the Hewlett Packard 6237B Triple output power 
supply is also set to DC 9V. The measurement is then commenced by getting the 
resonance frequency from the Hewlett Packard 3562A Dynamic signal analyzer. The 
microphone switch board is set to the microphone which is closest to the driver. Then 
the resonance frequency and the temperature of resonator was measured. The 
procedure to acquire resonance frequency has been mentioned in subsection 2 in 
this section. Getting quality factor also follows the same procedure previously 
mentioned. After finding resonance frequency, a data run is commenced by 
hooking up the Hewlett Packard 33 120 A function generator to the input amplifier for 
driver and setting the frequency and voltage. The driving frequency should be the 
resonance frequency obtained by previous procedure and the voltage is 1 Vpp. The 
acoustic amplitude is measured with each microphone. With the driver still running, the 
switch board selects a particular microphone and the value of amplitude at each 
position is measured with the SR 530 Lock - In amplifier. At this time, if there is 
variation in the resonator temperature, it is noted. All data recorded from the 16 
microphones are scaled by the correction factor obtained from microphone calibration. 
The corrected data go to computer program to plot mode shape. This concludes data 
acquisition of one mode. 
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Procedure for Wave Guide Characterization Collection 



1. Date: / / 97 ( ~ ) 

2. Ambient Condition 

Room Temper : °C ( K) / Speed of Sound = m/s 

Ambient Pressure: mbar / Initial Temp of resonator: °C 

3. Configuration (Constriction) 



# of Const. : EA 


/ 


Size : 


Area Ratio : 






Position of Const. : 


° (CCW/CW) 


from driver 


Position of Driver : 


o 




4. Resonance Freq: 


/ Q: 


/ Mode: 



5. Source information 

Input : Vpp / Type of source : 



6. Data : WG 


with 
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Receiver 
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Corerected 


Temp of 
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Figure 3-9. Data collection form. 
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IV. RESULTS AND DISCUSSION 



The purpose of this chapter is to show the comparison between theory and 
experimental results. The result of resonance frequency measurements are discussed first 
and then the mode shapes. 

A. RESONANCE FREQUENCIES 



The measured resonance frequencies are determined from pole-zero analysis of the 
frequency response. Because the resonance frequencies were measured at different 
temperatures, they were corrected to 20°C in order to make comparison of data 
recorded at various temperatures with the theoretical resonance frequencies more 
straightforward. The equation for conversion is 



f 20 ~ f , 



measured . 



1 + - 



20 

273 



1 + 



1 measured 

273 



(4.1) 



The theoretical resonance frequencies were calculated for two different cases. One case 
was without end corrections and the other was with end corrections. All the results are 
listed in Tables 4-1 and 4-2. Tables 4-la and 4.2a list the resonance frequencies for the 
low (L) and high (H) modes for the various constriction lengths and area ratios. Tables 4- 
lb and 4-2b lists the temperatures at which the frequencies were measured. Tables 4-1 c 
and 4-2c list the 20°C-corrected frequencies. Tables a, b and c are identical for Tables 4-1 
and 4-2. Tables 4-1 d and 4-2d give the calculated resonance frequencies. Table 4-ld 
ignores end corrections. Table 4-2d includes end corrections as calculated from higher 
order non propagating modes. Tables 4-le and 4-2e list the percent error between 
measured and calculated resonance frequencies. 

It is evident that the inclusion of end corrections has little effect on the prediction 
for the high modes. In both cases the average error is approximately 0.15%. However, 
end corrections are very important for the low mode. This point is evident from Fig. 4-1 
where it is seen that inclusion of end corrections is important for the low mode only. 
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Figure. 4-1 shows the 20°C-corrected resonance frequencies as well the predicted results 
(without and with end corrections) as functions of area ratio. 

In this section it was shown that end corrections are an important factor to include 
when calculating resonance frequencies, especially for the low mode. In the next section, 
it is shown that they are also important when calculating mode shapes. 
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1. Comparison of Measured Resonance Frequency with Calculated Results 
without End Corrections 



AR 


0.1 


0.2 


0.3 


0.5 


0.7 


1.0 


size 


L 


H 


L 


H 


L 


H 


L 


H 


L 


H 




L/16 


308.2 


461.5 


347.8 


458.1 


373.5 


455.4 


405.4 


449.3 


421.9 


444.0 


435.6 


L/32 


328.7 


449.6 


368.7 


447.8 


391.7 


446.4 


416.5 


443.6 


426.9 


439.3 




L/64 


348.6 


443.1 


384.2 


441.4 


402.6 


440.7 


422.1 


439.3 


431.1 


438.6 





Table 4- la. Measured resonance frequencies for different constrictions. 



AR 


0.1 


0.2 


0.3 


0.5 


0.7 


1.0 


size 


L 


H 


L 


H 


L 


H 


L 


H 


L 


H 




L/16 


22.5 


22.5 


21.4 


21.4 


21.8 


21.8 


22.0 


22.0 


22.4 


22.4 


22.1 


L/32 


23.9 


23.9 


23.3 


23.3 


23.6 


23.6 


23.7 


23.7 


21.6 


21.6 




L/64 


24.0 


24.0 


22.3 


22.3 


22.3 


22.3 


22.3 


22.3 


23.5 


23.5 





Table 4-lb. Temperatures for each measurement from table 4-la. 



AR 


0.1 


0.2 


0.3 


0.5 


0.7 


1.0 


size 


L 


H 


L 


H 


L 


H 


L 


H 


L 


H 




L/16 


306.9 


459.5 


347.0 


457.0 


372.4 


454.0 


404.0 


447.8 


420.2 


442.2 


434.0 


L/32 


326.5 


446.6 


366.6 


445.3 


389.3 


443.7 


413.9 


440.8 


425.7 


438.1 




L/64 


346.2 


440.1 


382.7 


439.7 


401.0 


439.0 


420.5 


437.6 


428.5 


436.0 





Table 4-lc. Measured resonance frequencies corrected to 20°C. 



AR 


0.1 


0.2 


0.3 


0.5 


0.7 


1.0 


size 


L 


H 


L 


H 


L 


H 


L 


H 


L 


H 




L/16 


317.6 


459.2 


359.7 


456.1 


383.2 


453.1 


409.2 


447.2 


421.5 


442.0 


434.2 


L/32 


349.3 


445.9 


387.3 


444.5 


404.2 


443.1 


420.2 


440.2 


427.1 


437.7 




L/64 


381.1 


439.5 


407.7 


438.8 


417.8 


438.1 


426.5 


436.7 


430.1 


435.5 





Table 4- Id. Calculated resonance frequencies at 20°C, without end corrections. 



AR 


0.1 


0.2 


0.3 


0.5 


0.7 


1.0 


size 


L 


H 


L 


H 


L 


H 


L 


H 


L 


H 




L/16 


3.48 


0.07 


3.71 


0.18 


2.91 


0.21 


1.29 


0.13 


0.32 


0.04 


0.05 


L/32 


6.98 


0.16 


5.65 


0.19 


3.83 


0.14 


1.53 


0.13 


0.33 


0.10 




L/64 


10.1 


0.13 


6.54 


0.20 


4.19 


0.20 


1.44 


0.20 


0.38 


0.11 





Table 4-le. Percent error between measured and calcu 



ated resonance frequency. 
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2. Comparison of Measured Resonance Frequency with Calculated Results with 

End Corrections 



AR 


0.1 


0.2 


0.3 


0.5 


0.7 


1.0 


size 


L 


H 


L 


H 


L 


H 


L 


H 


L 


H 




L/16 


308.2 


461.5 


347.8 


458.1 


373.5 


455.4 


405.4 


449.3 


421.9 


444.0 


435.6 


L/32 


328.7 


449.6 


368.7 


447.8 


391.7 


446.4 


416.5 


443.6 


426.9 


439.3 




L/64 


348.6 


443.1 


384.2 


441.4 


402.6 


440.7 


422.1 


439.3 


431.1 


438.6 





Table 4-2a. Measured resonance frequencies for different constrictions. 



AR 


0.1 


0.2 


0.3 


0.5 


0.7 


1.0 


size 


L 


H 


L 


H 


L 


H 


L 


H 


L 


H 




L/16 


22.5 


22.5 


21.4 


21.4 


21.8 


21.8 


22.0 


22.0 


22.4 


22.4 


22.1 


L/32 


23.9 


23.9 


23.3 


23.3 


23.6 


23.6 


23.7 


23.7 


21.6 


21.6 




L/64 


24.0 


24.0 


22.3 


22.3 


22.3 


22.3 


22.3 


22.3 


23.5 


23.5 





Table 4-2b. Temperatures for each measurement from table 4-2a. 



AR 


0.1 


0.2 


0.3 


0.5 


0.7 


1.0 


size 


L 


H 


L 


H 


L 


H 


L 


H 


L 


H 




L/16 


306.9 


459.5 


347.0 


457.0 


372.4 


454.0 


404.0 


447.8 


420.2 


442.2 


434.0 


L/32 


326.5 


446.6 


366.6 


445.3 


389.3 


443.7 


413.9 


440.8 


425.7 


438.1 




L/64 


346.2 


440.1 


382.7 


439.7 


401.0 


439.0 


420.5 


437.6 


428.5 


436.0 





Table 4-2c. Measured resonance frequencies corrected to 20°C. 



AR 


0.1 


0.2 


0.3 


0.5 


0.7 


1.0 


size 


L 


H 


L 


H 


L 


H 


L 


H 


L 


H 




L/16 


305.9 


459.1 


346.0 


456.0 


371.7 


453.0 


403.8 


447.1 


419.5 


442.0 


434.2 


L/32 


326.2 


445.9 


366.2 


444.4 


388.7 


443.0 


414.1 


440.2 


424.9 


437.7 




L/64 


345.6 


439.5 


381.4 


438.8 


400.0 


438.1 


420.0 


436.8 


427.8 


435.5 




Table 4-2d. Calculated 


resonance frequencies at 20°C with end corrections. 


AR 


0.1 


0.2 


0.3 


0.5 


0.7 


1.0 


size 


L 


H 


L 


H 


L 


H 


L 


H 


L 


H 




L/16 


0.34 


0.08 


0.28 


0.21 


0.20 


0.21 


0.04 


0.15 


0.16 


0.05 


0.05 


L/32 


0.10 


0.16 


0.11 


0.19 


0.15 


0.15 


0.04 


0.14 


0.18 


0.09 




L/64 


0.17 


0.13 


0.34 


0.20 


0.24 


0.20 


0.12 


0.19 


0.15 


0.12 




Ta 


Me 4-2e. Percent error between measured anc 


calculated resonance : 


requency. 
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measured & calculated resonance freq vs area ratio 




Figure 4- la. Comparison of calculated and measured resonance frequencies versus area 
ratio for the L/16 constriction. 

measured & calculated resonance freq vs area ratio 




Figure 4-lb. Comparison of calculated and measured resonance frequencies versus area 
ratio for the L/32 constriction. 
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measured & calculated resonance freq vs area ratio 




Figure 4-lc. Comparison of calculated and measured resonance frequencies versus area 
ratio for the L/64 constriction. 

B. MODE SHAPES 

Figures 4-2a and b show the comparison of the measured and computed mode 
shapes for the L/16, 0.1 area ratio constriction. The measurements have been normalized 
to the highest value of pressure amplitude. The amplitude of the mode shape is 
determined from least squares fit to the measured data. The driver is located 45° from the 
constriction. Figure 4-2a shows the results for the high mode. It is seen that there is 
good agreement between the measured and predicted mode shapes, even though end 
corrections are not included in the calculation. The driver is located at a place where both 
modes should be excited. However, the modes are well separated with this constriction, 
so there is little excitation of the low mode. Single mode excitation is indicated by a high 
standing wave ratio from the mode shape. 

Figure 4-2b shows the mode shape for the low mode. The agreement between the 
measured and calculated mode shape is significantly worse than with the high mode. 
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Figures 4-2c and d show the mode shapes again but this time end corrections are 
taken into account. Comparing Fig. 4-2c and 4-2a, it is seen that end corrections have 
little effect on the high mode. However, they are very important for the low mode, just 
as was discussed in the previous section. Figure 4-2d shows that when end corrections 
are taken into account, the agreement between calculated and measured mode shape 
improves significantly. 

Figures 4-3 a and b show the mode shapes for the case when the driver is located 
45° from the L/32, 0.7 area ratio constriction. In the case the constriction is relatively 
short and the area ratio is relatively large. As a result, the modes are not well separated. 
So although the driver is excited at the resonance frequency of one mode, there is 
significant excitation of the other mode, as evidenced by the low standing wave ratio. 

Figure 4-3 a shows the mode shape when the resonator is driven at the high mode 
resonance frequency. End correction are ignored in the calculation. The agreement 
between measured and calculated mode shape is not as good as in the previous, well 
separated-mode case (Fig. 4-2). Figure 4-3b shows the results for the low mode. Here 
the agreement between theory and experiment is worse. In both cases, the standing wave 
ratio is noticeably lower that in the well-separated case. 

Figure 4-3c and d show the effect of including end corrections. There is some 
improvement in the agreement for the high mode case. There is considerable improvement 
in the agreement for the low mode, and some improvement for the high mode. Many other 
comparisons are provided in the Appendix. The cases presented in this chapter are 
representative. The overall agreement in mode shape is good even in the overlapping- 
mode case. 

The results of this chapter indicate that end corrections are important both for 
calculation of resonance frequency and mode shapes. Also end corrections are more 
important for the low mode. This is because the acoustic impedance of the low mode is 
low in comparison to that of the high mode. The addition of the effective impedance of 
the end corrections is of more consequence, then, to the low mode. 
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Mode Shape for Length L/16, Area Ratio 0.1, and f=461.5 




Figure 4-2a. Mode shape for the high mode when the driver is located 45° from the 
constriction having length of L/16 and area ratio of 0.1. End corrections are not included 
in the calculation. 



Mode Shape for Length L/16, Area Ratio 0.1, and f=308.2 




Figure 4-2b. Mode shape for the low mode when the driver is located 45° from the 
constriction having length of L/16 and area ratio of 0.1. End corrections are not included 
in the calculation. 
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Mode Shape for Length L/16, Area Ratio 0.1, and f=461.5 




Figure 4-2c. Mode shape for the high mode when the driver is located 45° from the 
constriction having length of L/16 and area ratio of 0.1. End corrections are included in 
the calculation. 



Mode Shape for Length L/16, Area Ratio 0.1, and f=308.2 




Figure 4-2d. Mode shape for the low mode when the driver is located 45° from the 
constriction having length of L/16 and area ratio of 0. 1 . End corrections are included. 
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Mode Shape for Length L/32, Area Ratio 0.7, and f=439.3 




Figure 4-3a. Mode shape for the high mode when the driver is located 45° from the 
constriction having length of L/32 and area ratio of 0.7. End corrections are not included. 



Mode Shape for Length L/32, Area Ratio 0.7, and f=426.9 




Figure 4-3b. Mode shape for the low mode when the driver is located 45° from the 
constriction having length of L/32 and area ratio of 0.7. End corrections are not included. 
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Mode Shape for Length L/32, Area Ratio 0.7, and f=439.3 




Figure 4-3c. Mode shape for the high mode when the driver is located 45° from the 
constriction having length of L/32 and area ratio of 0.7. End corrections are included 



Mode Shape for Length L/32, Area Ratio 0.7, and f=426.9 




Figure 4-3 d. Mode shape for the low mode when the driver is located 45° from the 
constriction having length of L/32 and area ratio of 0.7. End corrections are included 
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V. 



SUMMARY AND CONCLUSIONS 



The purpose of this thesis is to investigate frequency splitting and mode shapes in a 
constricted annular resonator. The measurements are intended to serve as preliminary 
verification of a model predicting the performance of an annular thermoacoustic prime 
mover. 

The constriction splits the degeneracy of the plane wave-like modes of a uniform 
cross section annular resonator. One mode has a pressure antinode at the constriction, the 
other a pressure node. The resonance frequency of the former is higher than that of the 
latter. The difference in resonance frequency depends upon the length and porosity of the 
constriction. The effects of dissipation cause the frequency responses of the individual 
modes to overlap. 

An equivalent plane wave model is developed to predict the resonance frequencies 
and mode shapes. The experimental resonance frequencies are determined from pole-zero 
analysis of the measured frequency response. 

It is shown that it is critical to take into account end corrections at the constriction 
to get good agreement between the measured and predicted results. End corrections are 
more important for the low mode than the high mode. 

In general the percent error between measured and predicted resonance 
frequencies is less than 0.2 % for high mode and 0.3 % for the low mode. The measured 
mode shape are also in good agreement with predictions. 

The results of this thesis indicate that the simple constricted annular resonator is 
well understood. The sensitivity of the resonance frequencies and mode shapes to end 
corrections indicates that the annular resonator is a good geometry for investigating these 
effects further. For instance, how does one model the end corrections for a thermoacoustic 
stack? Are the end corrections for a stack and a constriction having the same porosity 
substantially different? These questions are left for future work. 
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APPENDIX DATA TABLE & MODE SHAPE 



Experiment 1. 



Constriction length: L/ 16 , area ratio: 0.1 , angle from driver: 180° 


Resonance, frq.: 463.7 Hz 


Q.109 Mode: high Temp: 24.6°C 


Microphone 

j (e) 


Vrms 

(mV) 


Corrected 

Data 


Microphone 
# (0) 


Vrms 

(mV) 


Corrected 

Data 


1 (000) 


70.5 


70.5 


9 (120) 


13.1 


16.6 


2 (006) 


67.3 


69.7 


10 (150) 


33.6 


50.1 


3 (012) 


55.9 


71.5 


11 (180) 


79.4 


67.0 


4 (018) 


52.8 


70.5 


12 (210) 


45.4 


65.1 


5 (024) 


44.3 


69.4 


13 (240) 


33.1 


43.0 


6 (030) 


58.9 


69.1 


14 (270) 


5.8 


7.7 


7 (060) 


39.2 


53.1 


15 (300) 


27.4 


28.3 


8 (090) 


15.9 


20.8 


16 (330) 


48.6 


57.0 



Mode Shape for Length L/16, Area Ratio 0.1, and f=463.7 
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Experiment 2. 



Constriction length: ' 


J 16 , area ratio: 0.1 , angle from driver:90° 


Resonance, frq.: 308.8 Hz 


Q: 68.5 Mode: 1 


ow Temp: 24.5°C 


Microphone 
# @) 


Vrms 

(mV) 


Corrected 

Data 


Microphone 
# (9) 


Vrms 

(mV) 


Corrected 

Data 


1 (000) 


39.5 


39.5 


9 (120) 


35.6 


44.7 


2 (006) 


17.2 


17.6 


10 (150) 


19.8 


29.5 


3 (012) 


3.4 


4.3 


11 (180) 


9.5 


8.1 


4 (018) 


20.0 


26.0 


12 (210) 


9.7 


13.8 


5 (024) 


28.6 


44.8 


13 (240) 


26.2 


34.1 


6 (030) 


45.4 


53.0 


14 (270) 


35.6 


47.9 


7 (060) 


44.9 


60.6 


15 (300) 


54.9 


57.8 


8 (090) 


43.4 


57.3 


16 (330) 


49.7 


58.7 



Mode Shape for Length L/16, Area Ratio 0.1 , and f=308.8 
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Experiment 3. 



Constriction length: L/ 16 , area ratio: 0.1 , angle from driver:45° 


Resonance frq.: 461.5 Hz 


Q: 98 Mode: high Temp: 22.5°C 


Microphone 
# (9) 


Vrms 

(mV) 


Corrected 

Data 


Microphone 
# 0) 


Vrms 

(mV) 


Corrected 

Data 


1 (000) 


49.8 


49.8 


9 (120) 


9.4 


11.9 


2 (006) 


48.7 


50.4 


10 (150) 


24.2 


36.1 


3 (012) 


39.8 


50.9 


11 (180) 


57.4 


48.5 


4 (018) 


37.8 


50.5 


12 (210) 


32.6 


46.8 


5 (024) 


31.4 


49.2 


13 (240) 


23.7 


30.8 


6 (030) 


41.7 


48.9 


14 (270) 


4.1 


5.5 


7 (060) 


27.8 


35.5 


15 (300) 


20.0 


20.7 


8 (090) 


11.3 


14.8 


16 (330) 


35.3 


41.4 



Mode Shape for Length L/16, Area Ratio 0.1, and f=461.5 
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Experiment 4. 



Constriction length: ] 


L/16 , area ratio: 0.1 , angle from driver :45° 


Resonance frq.: 308.2 Hz 


Q: 58 Mode: 1 


ow Temp: 22.5°C 


Microphone 

i 0) 


Vrms 

(mV) 


Corrected 

Data 


Microphone 

i (e) 


Vrms 

(mV) 


Corrected 

Data 


1 (000) 


36.4 


36.4 


9 (120) 


33.0 


41.5 


2 (006) 


17.0 


17.4 


10 (150) 


18.6 


27.7 


3 (012) 


2.3 


2.9 


11 (180) 


9.2 


7.8 


4 (018) 


17.7 


23.0 


12 (210) 


8.9 


12.7 


5 (024) 


25.8 


40.4 


13 (240) 


24.2 


31.5 


6 (030) 


41.3 


48.2 


14 (270) 


33.0 


44.4 


7 (060) 


41.3 


55.7 


15 (300) 


50.8 


53.5 


8 (090) 


40.3 


53.2 


16 (330) 


46.4 


54.8 



Mode Shape for Length L/16, Area Ratio 0.1 , and f=308.2 
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Experiment 5. 



Constriction length: L/32 , area ratio: 0.1 , angle from driver: 180° 


Resonance frq.: 448.9 Hz 


Q: 111 Mode: high Temp: 22.9°C 


Microphone 

# ( 9 ) 


Vrms 

(mV) 


Corrected 

Data 


Microphone 
# (0) 


Vrms 

(mV) 


Corrected 

Data 


1 (000) 


69.6 


69.6 


9 (120) 


20.6 


26.0 


2 (006) 


67.2 


69.5 


10 (150) 


37.9 


56.5 


3 (012) 


55.2 


70.5 


11 (180) 


82.8 


69.9 


4 (018) 


51.7 


68.9 


12 (210) 


44.0 


63.1 


5 (024) 


42.9 


67.1 


13 (240) 


29.7 


38.5 


6 (030) 


56.1 


65.7 


14 (270) 


2.5 


3.3 


7 (060) 


32.9 


44.5 


15 (300) 


30.7 


31.7 


8 (090) 


7.6 


9.9 


16 (330) 


50.2 


59.0 



Mode Shape for Length L/32, Area Ratio 0.1, and f=448.9 
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Experiment 6. 



Constriction length: ] 


L/32 , area ratio: 0.1 , an] 


E*le from driver:90° 


Resonance frq.: 328.7 Hz 


Q: 69.4 Mode: 1 


ow Temp: 22.3°C 


Microphone 
# (0) 


Vrms 

(mV) 


Corrected 

Data 


Microphone 

* ( 9 ) 


Vrms 

(mV) 


Corrected 

Data 


1 (000) 


20.3 


20.3 


9 (120) 


28.7 


35.9 


2 (006) 


3.8 


3.9 


10 (150) 


15.0 


22.2 


3 (012) 


21.2 


26.7 


11 (180) 


4.0 


3.4 


4 (018) 


27.9 


36.4 


12 (210) 


10.9 


15.5 


5 (024) 


26.1 


40.8 


13 (240) 


24.8 


32.2 


6 (030) 


37.4 


43.6 


14 (270) 


31.6 


42.3 


7 (060) 


36.7 


49.7 


15 (300) 


46.0 


48.3 


8 (090) 


35.6 


45.9 


16 (330) 


38.9 


45.8 



Mode Shape for Length L/32, Area Ratio 0.1 , and f=328.7 
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Experiment 7. 



Constriction 


length: 1 


J 32 , area ratio: 0.1 , angle from driver:45° 


Resonance frq.: 449.6 Hz 


Q: 90.3 


Mode: high Temp: 23.8°C 


Microphone 
# (0) 


Vrms 

(mV) 


Corrected 

Data 


Microphone 

# ( 9 ) 


Vrms 

(mV) 


Corrected 

Data 


1 (000) 


44.7 


44.7 


9 (120) 


13.3 


16.8 


2 (006) 


41.7 


43.1 


10 (150) 


24.4 


36.4 


3 (012) 


35.1 


44.8 


11 (180) 


52.5 


44.3 


4 (018) 


33.2 


44.3 


12 (210) 


28.3 


40.6 


5 (024) 


27.4 


42.9 


13 (240) 


19.2 


24.9 


6 (030) 


35.9 


42.0 


14 (270) 


1.7 


2.3 


7 (060) 


21.0 


28.4 


15 (300) 


19.9 


20.6 


8 (090) 


4.9 


6.4 


16 (330) 


32.5 


38.2 



Mode Shape for Length L/32, Area Ratio 0.1, and f=449.6 
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Experiment 8. 



Constriction length: ] 


U 32 , area ratio: 0.1 , an] 


gle from driver:45° 


Resonance frq.: 328.7 Hz 


Q: 56 Mode: 


ow Temp: 23.8°C 


Microphone 
# (0) __ 


Vrms 

(mV) 


Corrected 

Data 


Microphone 
# (0) 


Vrms 

(mV) 


Corrected 

Data 


1 (000) 


21.6 


21.6 


9 (120) 


29.0 


36.3 


2 (006) 


2.5 


2.6 


10 (150) 


15.3 


22.7 


3 (012) 


20.7 


26.0 


11 (180) 


4.3 


3.7 


4 (018) 


28.1 


36.6 


12 (210) 


11.0 


15.7 


5 (024) 


26.1 


40.8 


13 (240) 


24.9 


32.4 


6 (030) 


37.6 


43.9 


14 (270) 


31.8 


42.6 


7 (060) 


37.0 


50.1 


15 (300) 


46.5 


48.8 


8 (090) 


35.9 


46.3 


16 (330) 


39.5 


46.5 



Mode Shape for Length L/32, Area Ratio 0.1, and f=328.7 
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Experiment 9. 



Constriction 


length: L/64 , area ratio: 0. 1 , angle from driver: 180° 


Resonance frq.: 443.2 Hz 


Q: 116 


Mode: high Temp: 22.3°C 


Microphone 
# (9) 


Vrms 

(mV) 


Corrected 

Data 


Microphone 
# (0) 


Vrms 

(mV) 


Corrected 

Data 


1 (000) 


72.8 


72.8 


9 (120) 


24.9 


31.4 


2 (006) 


67.7 


69.9 


10 (150) 


41.4 


61.7 


3 (012) 


57.5 


73.3 


11 (180) 


85.9 


72.6 


4 (018) 


53.6 


71.5 


12 (210) 


45.1 


64.7 


5 (024) 


43.9 


68.7 


13 (240) 


29.7 


38.5 


6 (030) 


56.9 


66.6 


14 (270) 


1.4 


1.9 


7 (060) 


31.1 


42.0 


15 (300) 


33.0 


34.1 


8 (090) 


4.0 


5.2 


16 (330) 


52.9 


62.2 



Mode Shape for Length L/64, Area Ratio 0.1 , and f=443.2 
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Experiment 10. 



Constriction 


length: ] 


U 64 , area ratio: 0.1 , an; 


gle from driver:90° 


Resonance frq.: 350.2 Hz 


Q: 41 Mode: 


ow Temp: 24.0°C 


Microphone 
# 0) 


Vrms 

(mV) 


Corrected 

Data 


Microphone 

£ ( 9 ) 


Vrms 

(mV) 


Corrected 

Data 


1 (000) 


7.1 


7.1 


9 (120) 


17.9 


22.4 


2 (006) 


8.5 


8.7 


10 (150) 


9.2 


13.6 


3 (012) 


13.8 


17.4 


11 (180) 


2.0 


1.7 


4 (018) 


15.8 


20.7 


12 (210) 


7.3 


10.4 


5 (024) 


14.6 


22.8 


13 (240) 


16.0 


20.8 


6 (030) 


21.2 


24.8 


14 (270) 


19.8 


26.4 


7 (060) 


21.9 


29.5 


15 (300) 


28.0 


29.3 


8 (090) 


21.8 


28.2 


16 (330) 


22.3 


26.2 



Mode Shape for Length L/64, Area Ratio 0.1, and f=350.2 
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Experiment 11. 



Constriction length: 1 


L/64 , area ratio: 0.1 , angle from driver:45° 


Resonance frq.: 443.1 Hz 


Q: 87 Mode: high Temp: 24.0°C 


Microphone 
# (0) 


Vrms 

(mV) 


Corrected 

Data 


Microphone 

# ( 9 ) 


Vrms 

(mV) 


Corrected 

Data 


1 (000) 


39.8 


39.8 


9 (120) 


13.9 


17.5 


2 (006) 


37.0 


38.2 


10 (150) 


22.7 


33.8 


3 (012) 


31.1 


39.7 


11 (180) 


47.0 


39.7 


4 (018) 


29.1 


38.8 


12 (210) 


24.7 


35.4 


5 (024) 


23.9 


37.4 


13 (240) 


16.1 


20.9 


6 (030) 


31.0 


36.3. 


14 (270) 


0.8 


1.1 


7 (060) 


16.8 


22.7 


15 (300) 


18.6 


19.2 


8 (090) 


2.2 


2.9 


16 (330) 


292 


34.3 



Mode Shape for Length L/64, Area Ratio 0.1 , and f=443.1 
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Experiment 12. 



Constriction length: ] 


U 64 , area ratio: 0.1 , an 


gle from driver:45° 


Resonance frq.: 348.6 Hz 


Q: 57 Mode: 1 


ow Temp: 24.0°C 


Microphone 
# ( 0 ) 


Vims 

(mV) 


Corrected 

Data 


Microphone 
# (0) 


Vrms 

(mV) 


Corrected 

Data 


1 (000) 


10.0 


10.0 


9 (120) 


25.2 


31.6 


2 (006) 


12.7 


13.0 


10 (150) 


12.8 


19.0 


3 (012) 


20.5 


25.9 


11 (180) 


2.1 


1.8 


4 (018) 


23.3 


30.5 


12 (210) 


11.1 


15.8 


5 (024) 


21.4 


33.5 


13 (240) 


23.4 


30.4 


6 (030) 


31.0 


36.3 


14 (270) 


28.6 


38.2 


7 (060) 


31.6 


42.5 


15 (300) 


39.9 


41.8 


8 (090) 


31.2 


40.4 


16 (330) 


31.5 


37.0 



Mode Shape for Length L/64, Area Ratio 0.1, and f=348.6 




68 





Experiment 13. 



Constriction 


length: L/ 16 , area ratio: 0.7 , angle from driver: 180° 


Resonance frq.: 444.2 Hz 


Q: 115 


Mode: high Temp: 21.9°C 


Microphone 
# (6) 


Vrms 

(mV) 


Corrected 

Data 


Microphone 
# (6) 


Vrms 

(mV) 


Corrected 

Data 


1 (000) 


69.3 


69.3 


9 (120) 


15.4 


19.4 


2 (006) 


67.6 


69.9 


10 (150) 


34.7 


51.7 


3 (012) 


55.9 


71.3 


11 (180) 


80.0 


67.6 


4 (018) 


53.1 


70.8 


12 (210) 


46.6 


66.9 


5 (024) 


44.5 


69.6 


13 (240) 


35.4 


45.9 


6 (030) 


58.7 


68.7 


14 (270) 


9.2 


12.2 


7 (060) 


36.8 


49.8 


15 (300) 


22.6 


23.3 


8 (090) 


13.0 


17.0 


16 (330) 


45 6 


53.6 



Mode Shape for Length L/16, Area Ratio 0.7, and f=444.2 
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Experiment 14. 



Constriction length: ] 


L/16 , area ratio: 0.7 , angle from driver: 90° 


Resonance frq.: 421.2 Hz 


Q: 113 Mode:! 


ow Temp:22.1°C 


Microphone 
# (0) 


Vims 

(mV) 


Corrected 

Data 


Microphone 
i (9) 


Vrms 

(mV) 


Corrected 

Data 


1 (000) 


18.3 


18.3 


9 (120) 


54.3 


68.2 


2 (006) 


8.2 


8.5 


10 (150) 


32.5 


48.4 


3 (012) 


1.7 


2.2 


11 (180) 


16.2 


13.7 


4 (018) 


9.2 


12.2 


12 (210) 


16.6 


23.7 


5 (024) 


14.0 


21.9 


13 (240) 


42.7 


55.5 


6 (030) 


26.2 


30.7 


14 (270) 


52.6 


69.7 


7 (060) 


45.0 


60.7 


15 (300) 


68.2 


70.6 


8 (090) 


56.8 


74.0 


16 (330) 


45.0 


52.8 



Mode Shape for Length L/16, Area Ratio 0.7, and f=421.2 
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Experiment 15. 



Constriction 


length: L/16 , area ratio: 0.7 , an 


gle from driver:45° 


Resonance frq.: 444.0 Hz 


Q: 93 


Mode: high Temp: 22.4°C 


Microphone 

* ( 3 ) 


Vrms 

(mV) 


Corrected 

Data 


Microphone 

i (9) 


Vrms 

(mV) 


Corrected 

Data 


1 


(000) 


41.6 


41.6 


9 (120) 


10.7 


13.5 


2 


(006) 


40.5 


41.9 


10 (150) 


21.5 


32.0 


3 


(012) 


33.1 


42.2 


11 (180) 


48.3 


40.8 


4 


(018) 


31.5 


42.0 


12 (210) 


27.6 


39.6 


5 


(024) 


26.3 


41.2 


13 (240) 


20.6 


26.7 


6 


(030) 


34.5 


40.4 


14 (270) 


5.4 


7.2 


7 


(060) 


21.4 


29.0 


15 (300) 


15.1 


15.6 


8 


(090) 


7.5 


9.8 


16 (330) 


28.2 


33.1 



Mode Shape for Length L/16, Area Ratio 0.7, and f=444.0 
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Experiment 16. 



Constriction 


length: L / 16 , area ratio: 0.7 , an 


gle from driver:45° 


Resonance frq.: 421.9 Hz 


Q: 117 


Mode: low Temp: 22.4°C 


Microphone 


Vrms 


Corrected 


Microphone 


Vrms 


Corrected 


# (0) 


(mV) 


Data 


i (e) 


(mV) 


Data 


1 (000) 


14.3 


14.3 


9 (120) 


43.4 


54.5 


2 (006) 


6.7 


6.9 


10 (150) 


25.8 


38.4 


3 (012) 


3.1 


3.9 


11 (180) 


12.7 


10.8 


4 (018) 


8.4 


11.2 


12 (210) 


14.4 


20.6 


5 (024) 


11.9 


18.6 


13 (240) 


34.9 


45.3 


6 (030) 


21.9 


25.6 


14 (270) 


42.4 


56.2 


7 (060) 


36.7 


49.5 


15 (300) 


54.6 


56.5 


8 (090) 


45.9 


59.8 


16 (330) 


35.7 


41.9 



Mode Shape for Length L/16, Area Ratio 0.7, and f=421.9 
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Experiment 17. 



Constriction length: L/32 , area ratio: 0.7 , angle from driver: 180° 


Resonance frq.: 441.1 Hz 


Q: 114 Mode: high Temp: 23.8°C 


Microphone 
# (9) 


Vrms 

(mV) 


Corrected 

Data 


Microphone 
# (9) 


Vrms 

(mV) 


Corrected 

Data 


1 (000) 


69.4 


69.4 


9 (120) 


20.4 


25.7 


2 (006) 


67.4 


69.6 


10 (150) 


37.5 


55.9 


3 (012) 


55.5 


70.8 


11 (180) 


81.4 


68.8 


4 (018) 


51.5 


68.7 


12 (210) 


45.5 


65.2 


5 (024) 


43.6 


68.2 


13 (240) 


32.4 


42.0 


6 (030) 


56.8 


66.5 


14 (270) 


5.7 


7.6 


7 (060) 


32.8 


44.3 


15 (300) 


26.7 


27.6 


8 (090) 


7.8 


10.2 


16 (330) 


48.0 


56.4 



Mode Shape for Length L/32, Area Ratio 0.7, and f=441.1 
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Experiment 18. 



Constriction 


length: ] 


U 32 , area ratio: 0.7 , angle from driver:90° 


Resonance frq.: 426.2 Hz 


Q: 129 


Mode: 1 


ow Temp:21.1°C 


Microphone 


Vrms 


Corrected 


Microphone 


Vrms 


Corrected 


# (0) 


(mV) 


Data 


# (0) 


(mV) 


Data 


1 (000) 


8.7 


8.7 


9 (120) 


57.8 


72.6 


2 (006) 


2.8 


2.9 


10 (150) 


31.7 


47.2 


3 (012) 


10.7 


13.6 


11 (180) 


7.2 


6.1 


4 (018) 


17.2 


22.9 


12 (210) 


24.6 


35.2 


5 (024) 


20.3 


31.8 


13 (240) 


52.1 


67.6 


6 (030) 


34.3 


40.1 


14 (270) 


59.4 


78.8 


7 (060) 


52.8 


71.3 


15 (300) 


71.9 


74.4 


8 (090) 


63.5 


82.8 


16 (330) 


42.4 


49.8 



Mode Shape for Length L/32, Area Ratio 0.7, and f=426.2 
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Experiment 19. 



Constriction length: ] 


U 32 , area ratio: 0.7 , angle from driver:45° 


Resonance frq.: 439.3 Hz 


Q: 101 Mode: high Temp: 21.6°C 


Microphone 
# (0) 


Vims 

(mV) 


Corrected 

Data 


Microphone 

# (B) 


Vrms 

(mV) 


Corrected 

Data 


1 (000) 


43.4 


43.4 


9 (120) 


14.6 


18.4 


2 (006) 


42.0 


43.4 


10 (150) 


24.0 


35.8 


3 (012) 


34.3 


43.8 


11 (180) 


51.6 


43.6 


4 (018) 


31.7 


42.3 


12 (210) 


27.9 


40.0 


5 (024) 


26.7 


41.8 


13 (240) 


19.8 


25.7 


6 (030) 


34.8 


40.7 


14 (270) 


5.8 


7.7 


7 (060) 


20.2 


27.3 


15 (300) 


18.8 


19.4 


8 (090) 


7.0 


9.1 


16 (330) 


30.6 


36.0 



Mode Shape for Length L/32, Area Ratio 0.7, and f=439.3 
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Experiment 20. 



Constriction length: ] 


U 32 , area ratio: 0.7 , an 


gle from driver:45° 


Resonance frq.: 426.9 Hz 


Q: 125 Mode: 1 


ow Temp:21.6°C 


Microphone 
£ (6) 


Vrms 

(mV) 


Corrected 

Data 


Microphone 

i (B) 


Vrms 

(mV) 


Corrected 

Data 


1 (000) 


8.7 


8.7 


9 (120) 


42.3 


53.1 


2 (006) 


6.5 


6.7 


10 (150) 


23.2 


34.5 


3 (012) 


9.5 


12.1 


11 (180) 


9.6 


8.1 


4 (018) 


13.7 


18.2 


12 (210) 


19.1 


27.4 


5 (024) 


15.7 


24.5 


13 (240) 


38.7 


50.2 


6 (030) 


26.1 


30.5 


14 (270) 


43.6 


57.8 


7 (060) 


39.2 


52.9 


15 (300) 


52.5 


54.3 


8 (090) 


46.7 


60.9 


16 (330) 


30.9 


36.3 



Mode Shape for Length L/32, Area Ratio 0.7, and f=426.9 
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